Charge transport of electron doped Mott insulators on a triangular lattice 
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The charge transport of electron doped Mott insulators on a triangular lattice is investigated within 
the t-J model. The conductivity spectrum shows a low-energy peak and an unusual midinfrared 
band, while the resistivity is characterized by a crossover from the high temperature metallic-like 
to low temperature insulating-like behavior, in qualitative agreement with experiments. Our results 
also show that the mechanism that causes this unusual charge transport is closely related to a 
competition between the kinetic energy and magnetic energy in the system. 
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The recent discovery of superconductivity in doped 
cobaltates, Naa;Co02 • yH20 with the superconducting 
transition temperature Tc ^ 5K, has generated great in- 
terests due to the role of strong electron correlations^. 
This compound has a lamellar structure consisting of 
the two-dimensional (2D) C0O2 layers separated by a 
thick insulating layer of Na"*" ions and H2O molecules, 
where the one-half spin Co"*^ ions are arranged on a 
triangular lattice. This structure is similar to cuprates 
superconductors in the sense that they also have a lay- 
ered structure of the square lattice of the Cu02 plane 
separated by insulating layers^. It has been well estab- 
lished that Cu2+ ions exhibit an antiferromagnetic (AF) 
long-range order (AFLRO) in the parent compounds of 
cuprate superconductors, where superconductivity oc- 
curs when the AFLRO state is suppressed by hole or elec- 
tron doping^''^. However, Na2;Co02-yH20 is viewed as an 
electron doped Mott insulator, where superconductivity 
appears with electron doping*. A fundamental similar- 
ity between cuprate and cobaltate superconductors has 
been seen in the decreases in the superconducting tran- 
sition temperature for both underdoped and overdoped 
materials*. The optimal doping for superconductivity 
occurs at 0.3 electrons per Co above the ground-state 
NaoCo02 • I.3H2O, which is a half-filled two-electron t2g 
derived band, while for cuprates, the optimal doping oc- 
curs at 0.15 holes (electrons) added to the half-filled band 
of the parent compound^. On the other hand, doped 
Mott insulators on a triangular lattice, where the geo- 
metric spin frustration exists, are also of interests in their 
own right, with many unanswered fascinating questions. 
Historically the undoped Mott insulator on a triangular 
lattice was firstly proposed to be a model for microscopic 
realization of the spin liquid due to the existence of the 
strong spin frustration^. It has been argued that this spin 
liquid state may play a crucial role in the mechanism for 
superconductivity in doped cuprates^. Thus the unex- 
pected finding of superconductivity in doped cobaltates 
has raised the hope that it may help solve the unusual 
physics in doped cuprates. 

Some experimental studies have revealed a non-Fermi 
liquid behavior in doped cobaltates*'^'*. Among the 



striking features in the normal-state, a hallmark is the 
charge transport, where the resistivity shows a tempera- 
ture linear variation in a wide range of temperatures*'^. 
Moreover, this unusual charge transport is found to be 
intriguingly related to the magnetic correlation. The 
charge transport measurement is a powerful probe for 
interacting electron systems^, and can provide very de- 
tailed knowledges about the low-energy excitations as 
electrons are doped to cobaltates. One of the impor- 
tant issues in experimental and theoretical investigations 
of doped cobaltates is understanding of the low-energy 
excitations in these compounds, where the clue to their 
superconductivity lies possibly already in their noncon- 
ventional normal-state properties. Therefore it is of in- 
terest to have a detailed look at the charge transport. In 
this paper, we try to study this issue. We show that the 
resistivity of electron doped cobaltates is characterized 
by a crossover from the high temperature metallic-like to 
low temperature insulating-like behavior, and the mecha- 
nism that causes this unusual transport is closely related 
to a competition between the kinetic energy and mag- 
netic energy in the system. 

In electron doped cobaltates, the characteristic fea- 
ture is the presence of the 2D C0O2 plane* as mentioned 
above, and it seems evident that the unusual behaviors 
are dominated by this plane. It has been argued that the 
essential physics of the doped C0O2 plane is contained in 
the t-J model on a triangular lattice*". 
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where > 0, the summation is over all sites i, and for 
each i, over its nearest-neighbor 77, Cj^ (Cic) is the elec- 
tron creation (annihilation) operator, Sj — CjaCi/2 is 
the spin operator with cr = {ax, o'y,az) as the Pauli ma- 
trices, /X is the chemical potential, and the projection 
operator P removes zero occupancy, i.e., n^j + > 1 
with Uicr = Cj^Cia- In the past fifteen years, some useful 
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methods have been proposed to treat the no double oceu- 
pancy local constraint in hole doped euprates^^. In par- 
ticular, a fermion-spin theory based on the partial charge- 
spin separation has been developed to study the physical 
properties of doped cuprates^^, where the no double oc- 
cupancy local constraint can be treated properly in ana- 
lytical calculations. Within this theory, the charge trans- 
port of the underdoped cuprates has been discussed^^'^^, 
and the results are consistent with experiments'''^'*. To 
apply this theory in electron doped cobaltates, the t-J 
model (1) can be rewritten in terms of a particle- hole 
transformation Cjo- — > //_^ as 
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H = -tY. W + E + E Si • Si+^, (2) 



supplemented by a local constraint flafi^ < 1 to 
remove double occupancy, where as a matter of conve- 
nience, we have set t = —te, (fia) is the hole creation 
(annihilation) operator, and Sj = fjafi/2 is the spin 
operator in the hole representation. Now we follow the 
partial charge-spin separation fermion-spin theory^^, and 
decouple hole operators /j-f and fii as. 
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where the spinful fermion operator Uicr = e~**'<'ai 
describes the charge degree of freedom together with 
the phase part of the spin degree of freedom (dressed 
fermion), while the spin operator Si describes the am- 
plitude part of the spin degree of freedom (dressed 
spinon). then the no double occupancy local constraint, 

S~ Sf ) = 1 — a\ai < 1, is satisfied in analyti- 
cal calculations, and the double dressed fermion occu- 
e^^'^alale**'-" = and ai^ai-^ = 
= 0, are ruled out automatically. It 
has been shown that these dressed fermion and spinon 
are gauge invariant, and in this sense, they are real and 
can be interpreted as the physical excitations. Since the 
phase factor e"'**" is separated from the bare spinon 
operator, therefore it describes a spinon cloud^^. In this 
case, the dressed fermion Qji(j IS 3; spinless fermion in- 
corporated a spin cloud e"'*'" (magnetic fiux), and is a 
magnetic dressing. In other words, the gauge invariant 
dressed fermion carries some spin messages, i.e., it shares 
some effects of the spin configuration rearrangements due 
to the presence of the electron itself^^. Although in com- 
mon sense is not a real spinful fermion, it behaves like 
a spinful fermion. In this partial charge-spin separation 
fermion-spin representation, the low-energy behavior of 
the t-J model (2) can be expressed as^^, 
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with Jeff = (1 — (^)^J, and S = {a\^aic) = (alaj) is the 
electron doping concentration. At the zero doping, the 
t-J model is reduced to the Heisenberg model. Many 
authors^^ have shown unambiguously that as in a square 
lattice, there is indeed AFLRO in the ground state of the 
AF Heisenberg model on a triangular lattice. However, 
this AFLRO is destroyed more rapidly with increasing 
doping than that on a square lattice due to the strong 
geometry frustration. Therefore there is no AFLRO away 
from the zero doping, i.e., {8^) = 0. On the other hand, 
in the electron underdoped regime where superconductiv- 
ity occurs, the weak ferromagnetism can be induced*^' 
since the effective magnetic coupling in the present case 
is Jeff — 2(f)\t\ with the dressed fermion mean-field (MF) 
particle- hole parameter cf) = {al^ai^fjrr) ■ 

Since the no double occupancy local constraint has 
been treated properly within the partial charge-spin sep- 
aration fermion-spin theory, this leads to disappearing of 
the extra gauge degree of freedom related to the local 
constraint"'^^. In this case, the charge fluctuation couples 
only to dressed fermions^^'^^. Based on this theory, the 
charge transport of the hole doped Mott insulator on a 
square lattice has been discussed-'^^'-'^'', and the results are 
consistent with experiments^ Following their discus- 
sions, the optical conductivity of electron doped Mott 
insulators on a triangular lattice can be obtained as, 

a(c.) = {Zxtef^Y.^l,p ^4")(k,.;'+u;) 
x4")(k,a;0 "^('"^ + ^^-"^(^^ (5) 

UJ 

where Z is the number of the nearest neigh- 
bor sites, the spinon correlation function x — 
{8+Sr^f), 7?fe - {[sink, + sin(fc,/2)cos(V3fc,/2)]2 + 
3[sin(-\/3A;y/2)cos(A;x/2)]^}/9, npi^o) is the fermion dis- 
tribution function, and the dressed fermion spectral func- 
tion Ai^\k ,Lu) is obtained as A^a^\k,Lu) = —2lmg^{k,Lu). 
The full dressed fermion Green's function g~^[k,uj) = 
gi"^^^{k,uj) - 'E^°-'>{k,u>) with the MF dressed fermion 

Green's function 9^"'' ^ {k,u}) = to — ^k, and the second- 
order dressed fermion self-energy from the dressed spinon 
pair bubble^^'^^. 
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2 Z_/ 4iJg+pLVg 

F^'\k,p,q) 

UJ + ^q+p ~ — ^p+fe 

F^^\k,p,q) 

W — OJq+p + Ulq — ^p+k 
U! + LOq+p + LOq — £.p+k 

F^^\k,p,q) \ 

W — ^q+p — <^q— ^p+k J 
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where Bk = X[{2ex' +xhk- (ex + 2x')], A = 2ZJ,s, e = 
l + 2t(j)/JeB, 7/c = [coskx + 2 COS {kx/ 2) COS {V^ky/ 2)]/ 3, 
and 

F^^\k,p,q) =nF{£,p+k)[nB{uJq) - nB{uJq+p)\ 

+ nB{iOq+p)[l + ns(a;g)], (7a) 
F^'^\k,p,q) = nF(Cp+fe)[ni3K+p) - "sH)] 

+ nB{LJq)[l + nB{c^q+p)], (7b) 
F(3) (fc, p, q) = np (Cp+fc ) [1 + ns (w<,+p) + ns (wg)] 

+ nB(t^g)nB(t^g+p), (7c) 
F^^\k,p,q) = [l + nB{iOqm + nB{LOq+p)] 

- »^F(^p+fe)[l + nB{u)g+p) + nsiuiq)], (7d) 

with nB(u;) is the Bose distribution function, and the MF 
dressed fermion and spinon excitation spectra are given 
by, 



Ztxik - 

^17^ + ^27/0+ A3, 



respectively, where 

Ai=ae\^(^X + ^x"), 
A2 = eA^{l-Z)a{\ex + xl 
-a(C^ + ic)-^(l-a)], 



X^[a{C^ + le^C) + ^{l-a){l + e^) 
ae-^(^X + ex')], 



(8a) 
(8b) 
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with spinon correlation 

functions ^ = {SfSf^^, C = (1/^^) E^v(^i^^r+,,), 
and = {l/Z^)j:fi^,{S!+f,Sl^,). In order to satisfy 

the sum rule for the correlation function {S^S^) = 1/2 
in the absence of AFLRO, a decoupling parameter a has 
been introduced in the MF calculation, which can be re- 
garded as the vertex correction^^'^''. All these mean-field 
order parameters, decoupling parameter, and the chemi- 
cal potential are determined self-consistently^*. 

We have performed a numerical calculation for the op- 
tical conductivity in Eq. (5), and the results at elec- 
tron doping S = 0.40 (solid line), S = 0.35 (dashed line), 
and 6 = 0.30 (dotted line) for parameter t/J = —2.5 
in temperature T = are shown in Fig. 1, where the 
charge e has been set as the unit. The conductivity spc;c;- 
trum shows a low-energy peak appearing at a; ~ 0, which 
decays rapidly, and a broad midinfrared band. This 
midinfrared band is electron doping dependent, and the 
component increases with increasing electron doping for 
l|t|<a;<3|f|, and is nearly independent of elec- 
tron doping for (jj > 3 | t |. This reflects an increase in 
the mobile carrier density, and indicates that the spec- 
tral weight of the midinfrared sideband is taken from the 



Drude absorption, then the spectral weight from both low 
energy peak and midinfrared band represent the actual 
free-carrier density. For a better understanding of the 
optical properties of electron doped Mott insulators on 
a triangular lattice, we have studied the conductivity at 
different temperatures, and the results at electron doping 
5 = 0.35 for t/J = —2.5 with temperature T = (solid 
line), T = 0.3J (dashed hue), and T = 0.7 J (dotted line) 
are plotted in Fig. 2. It is shown that the conductivity 
spectrum is temperature-dependent for w < 3 | t |, and 
almost temperature-independent for a; > 3 | f | . The low- 
energy peak broadens and decreases in the height with 
increasing temperatures, and there is a tendency towards 
the Drude-like behavior. At the same time, the midin- 
frared band is severely suppressed with increasing tem- 
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FIG. 1. The conductivity of electron doped Mott insulators 
on a triangular lattice as a function of frequency at electron 
doping 5 = 0.40 (solid line), 5 — 0.35 (dashed line), and 
S = 0.30 (dotted line) for t/J = —2.5 with temperature T = 0. 
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FIG. 2. The conductivity of electron doped Mott insulators 
on a triangular lattice as a function of frequency at electron 
doping 5 = 0.35 for t/J= —2.5 with temperature T — (solid 
line), T = 0.3J (dashed line), and T = 0.7J (dotted line). 
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peratures, and vanishes at high temperatures, in quah- 
tative agreement with these of hole doped Mott insula- 
tors on a square lattice^'^^'"'^^'^^. As in doped cuprates, 
the charge transport is governed by the dressed fermion 
scattering, therefore 6 dressed fermion are responsible for 
the conductivity, i.e., the optical conductivity in elec- 
tron doped cobaltates is carried by 6 electrons. Since 
the strong electron correlation is common for both hole 
doped cuprates and electron doped cobaltates, these sim- 
ilar behaviors observed from the optical conductivity are 
expected. 

Now we tiirn to discuss the resistivity, which can 
be obtained in terms of the conductivity as p(T) = 
1/ lim^^o and the results for t/J = —2.5 at elec- 

tron doping S = 0.36 (solid line), S = 0.35 (dashed 
line), and S = 0.34 (dotted line) are plotted in Fig. 3, 
in comparison with the experimental data^ taken from 
NaxCo02 • yH2 (inset). Our present results show obvi- 
ously that p{T) is characterized by a crossover from the 
high temperature metallic-like {dp{T)/dT > 0) to low 
temperature insulating- like behavior {dp{T)/dT < 0), 
but the metallic-like linear temperature dependence dom- 
inates over a wide temperature range, which is qualita- 
tively consistent with experiments^'^. 

An explanation for this unusual charge transport can 
be found from a competition between the kinetic energy 
and magnetic energy, since the present t-J model (4) is 
characterized by the competition between the kinetic en- 
ergy {6t) and magnetic energy (J), with the magnetic 
energy J favors the magnetic order for spins, while the 
kinetic energy St favors dclocalization of electrons and 
tends to destroy the magnetic order. Although both 
dressed fermions and spinons contribute to the charge 
transport in the partial charge-spin separation fermion- 



12 




4 ' ' ' ' ' 

0.00 0.05 0.10 0.tS 0.20 

T/J 

FIG. 3. The resistivity of electron doped Mott insula- 
tors on a triangular lattice as a function of temperature for 
t/J = —2.5 at electron doping 5 = 0.36 (solid line), S = 0.35 
(dashed line), and 6 = 0.34 (dotted line). Inset: the experi- 
mental result on Nao.3BCo02 • I.3H2O taken from Ref. [1]. 



spin theory, the dressed fermion scattering dominates the 
charge transport^^'^'^. The dressed fermion scattering 
rate is obtained from the full dressed fermion Green's 
function (then the dressed fermion self-energy and spec- 
tral function) by considering the interaction between the 
dressed fermion and spinon. In this case, the crossover 
from the high temperature metallic-like to low tempera- 
ture insulating-like behavior in the resistivity of electron 
doped cobaltates is closely related to this competition. In 
lower temperatures, the dressed fermion kinetic energy is 
much smaller than the magnetic energy, then the mag- 
netic fluctuation is strong enough to severely reduce the 
dressed fermion scattering and thus is responsible for the 
insulating-like behavior in the resistivity. With increas- 
ing temperatures, the dressed fermion kinetic energy is 
increased, while the dressed spinon magnetic energy is 
decreased. In the region where the dressed fermion ki- 
netic energy is much larger than the dressed spinon mag- 
netic energy at higher temperatures, the dressed fermion 
scattering would give rise to the metallic-like resistivity. 

We emphasize that this competition between the ki- 
netic energy and magnetic energy exists in all doped Mott 
insulators. However, in the doped two-leg ladder Mott in- 
sulators, the charged carrier's motion is also suppressed 
by interference effects between the two legs^^, therefore 
in the doped two-leg ladder Mott insulators both com- 
petition between the kinetic energy and magnetic energy 
and interference effects between the two legs cause the 
unusual charge transport. This is why even in the hole 
doped two-leg ladder systems, the metallic-like state ap- 
pears in much higher doping concentration. On the other 
hand, hole doped Mott insulators (opposite sign of hop- 
ping t) on a triangular lattice have been discussed by 
many authors^^'^'', where the conventional quasiparticlc 
picture may be broken by the effects of geometric spin 
fluctuation^^. These and related differences^^ between 
hole and electron doped Mott insulators on a triangu- 
lar lattice reflect that there is no particle-hole symmetry, 
and the sign of t is important^". 

In summary, we have studied the charge transport of 
electron doped cobaltates within the t-J model based on 
the partial charge-spin separation fermion-spin theory. 
It is shown that the conductivity spectrum shows a low- 
energy peak and a broad midinfrared band, while the 
resistivity is characterized by a crossover from the high 
temperature metallic-like to low temperature insulating- 
like behavior. Our results also show that the mechanism 
that causes this unusual charge transport is closely re- 
lated to the competition between the kinetic energy and 
magnetic energy in the system. 

Recent experimental studies^'' have shown that the su- 
perconducting transition temperature in Naa,Co02 -2/1120 
is proportional to the electron doping concentration in 
the underdoped regime, and satisfies the Uemura rela- 
tion for hole doped cuprates. This and other exper- 
iments have shown that there is a remarkable resem- 
blance in superconducting-state properties between the 
electron doped cobaltate Naa;Co02-yH20 and hole doped 
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cupratcs 



,24.4,1 



manifesting that two systems have simi- 



lar underlying superconducting mechanism. Within the 
partial charge-spin separation fermion-spin theory, the 
mechanism of superconductivity in hole doped cuprates^^ 
has been discussed, where dressed holons interact occur- 
ring directly through the kinetic energy by exchanging 
the dressed spinon excitations, leading to a net attrac- 
tive force between the dressed holons, then the electron 
Cooper pairs originating from the dressed holon pairing 
state are due to the charge-spin recombination, and their 
condensation reveals the superconducting ground-state. 
In this case, the electron superconducting transition tem- 
perature is determined by the dressed holon pair transi- 
tion temperature, and is proportional to the hole doping 
concentration in the underdoped regime, in agreement 
with experiments. Since the strong electron correlation 
is common for both hole doped cupratcs and electron 
doped cobaltates, then it is possible that superconduc- 
tivity in electron doped cobaltates is also driven by the 
kinetic energy as in hole doped cupratcs. Following the 
discussions in Ref.^^, we have studied this issue^^, and 
found that the superconducting transition temperature 
in electron doped cobaltates is suppressed to a lower tem- 
perature due to the strong frustration. These and related 
theoretical results will be presented elsewhere. 

After submitting this paper, we became aware of re- 
cent optical measurements^'' on Nao.7Co02 supporting 
our theoretical results. 



ACKNOWLEDGMENTS 

The authors would like to thank Dr. Jihong Qin, Dr. 
Tianxing Ma, and Professor Y.J. Wang for the helpful 
discussions. This work was supported by the National 
Natural Science Foundation of China under Grant Nos. 
10125415, 90103024, and 10347102, and the National Sci- 
ence Council. 



* To whom correspondence should be addressed. 

^ K. Takada, H. Sakural, E. Takayama-Muromachi, F. Izuml, 
R.A. Dilanian, and T. Sasaki, Nature 422, 53 (2003). 

^ M.A. Kastner, R.J. Birgeneau, G. Shiran, and Y. Endoh, 
Rev. Mod. Phys. 70, 897 (1998). 

^Y. Tokura, H. Takagi, and S. Uchida, Nature 337, 345 
(1989); H. Takagi, S. Uchida, and Y. Tokura, Phys. Rev. 
Lett. 62, 1197 (1989). 

" R.E. Schaak, T. Klimczuk, M.L. Foo, and R.J. Cava, Na- 
ture 424, 527 (2003). 

^ P.W. Anderson, Mater. Res. Bull, 8, 153 (1973); P. Fazekas 
and P.W. Anderson, Philos. Mag. 30, 423 (1974); T.K. Lee 
and Shiping Feng, Phys. Rev. B 41, 11110 (1990). 

® P.W. Anderson, in Frontiers and Borderlines in Many Par- 



ticle Physics, edited by R.A. Broglia and J.R. Schrieffer 
(North-Holland, Amsterdam, 1987), p. 1; Science 235, 
1196 (1987). 

^ Y. Wang, N.S. Rogado, R.J. Cava, and N.P. Ong, Nature 

423, 425 (2003). 
® F.C. Chou, J.H. Cho, P.A. Loo, E.T. Abel, K. Matan, and 

Y.S. Lee, cond-mat/0306659, unpublished. 
® D.B. Tanner and T. Timusk, in Physical Properties of High 

Temperature Superconductors III, edited by D.M. Ginsberg 

(World Scientific, Singapore, 1992), p. 363. 
1° G. Baskaran, Phys. Rev. Lett. 91, 097003 (2003); B. Ku- 
mar and B.S. Shastry, Phys. Rev. B 68, 104508 (2003); Q. 

Wang, D.H. Lee, and P.A. Lee, Phys. Rev. B 69, 092504 

(2003). 

See, e.g., Proceedings of Los Alamos Symposium, edited 
by K.S. Bedell, D. Coffey D.E. Mehzer, D. Pines, and 
J.R. Schrieffer (Addison- Wesley, Redwood city, California, 
1990). 

Shiping Feng, Jihong Qin, and Tianxing Ma, J. Phys.: Con- 
dens. Matter 16, 343 (2004); Shiping Feng, Tianxing Ma, 
and Jihong Qin, Mod. Phys. Lett. B 17, 361 (2003); Ship- 
ing Feng, Z.B. Su, and L. Yu, Phys. Rev. B 49, 2368 (1994). 
Shiping Feng and Zhongbing Huang, Phys. Lett. A 232, 
293 (1997); Feng Yuan, Jihong Qin, Shiping Feng, and Wei 
Yen Chen, Phys. Rev. B 67, 134504 (2003). 

"'^ Y. Ando, A.N. Lavrov, S. Komiya, K. Segawa, and X.F. 
Sun, Phys. Rev. Lett. 87, 017001 (2001). 

1^ G.B. Martins, R. Eder, and E. Dagotto, Phys. Rev. B60, 
R3716 (1999); G.B. Martins, J.C. Xavier, C. Gazza, M. 
Vojta, and E. Dagotto, Phys. Rev. B63, 014414 (2000); 
G.B. Martins, C. Gazza, J.C. Xavier, A. Feiguin, and E. 
Dagotto, Phys. Rev. Lett. 84, 5844 (2000). 

^® L. Caprioti, A.E. Trumper, and S. Sorella, Phys. Rev. Lett. 
82, 3899 (1999). 

"P.W. Anderson, Science 288, 480 (2000); cond- 
mat/0108522, unpubUshed. 

Shiping Feng and Yun Song, Phys. Rev. B 55, 642 (1997). 
J. Kondo and K. Yamaji, Prog. Theor. Phys. 47, 807 
(1972). 

^° W. Stephan and P. Horsch, Phys. Rev. B 42, 8736 (1990); 
A. Moreo and E. Dagotto, Phys. Rev. B 42, 4786 (1990). 
Jihong Qin, Yun Song, Shiping Feng, and Wei Yeu Chen, 
Phys. Rev. B 65, 155117 (2002). 

L.O. Manuel, C.J. Gazza, A.E. Feiguin, and A.E. Trumper, 
J. Phys.: Condens. Matter 15, 2435 (2003); A.E. Trumper, 
C.J. Gazza, and L.O. Manuel, cond-mat/0306147, unpub- 
lished. 

W.Q. Yu, Shiping Feng, Z.B. Huang, and H.Q. Lin, Phys. 
Rev. B 61, 7409 (2000). 

A. Kanigel, A. Keren, L. Patlagan, K.B. Chashka, B. 
Fisher, P. King, and A. Amato, cond-mat/0311427, un- 
published. 

Shiping Feng, Phys. Rev. B 68, 184501 (2003). 
Bin Liu, Ying Liang, and Shiping Feng (unpublished). 
N.L. Wang, P. Zheng, D. Wu, Y.C. Ma, T. Xiang, R.Y. 
Jin, and D. Mondrus, cond-mat/0312630, unpublished. 



5 



